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ABSTRACT 

We study the variation of the frequency sphtting coefficients describing the solar as- 
phericity in both GONG and MDI data, and use these data to investigate temporal 
sound-speed variations as a function of both depth and latitude during the period from 
1995-2000 and a little beyond. The temporal variations in even splitting coefficients 
are found to be correlated to the corresponding component of magnetic flux at the 
solar surface. We confirm that the sound-speed variations associated with the surface 
magnetic field are superficial. Temporally averaged results show a significant excess in 
sound speed around r — 0.92Rq and latitude of 60°. 

Key words: Sun:Interior - Sun:Oscillations 



6 



1 INTRODUCTION 

Helioseismology - the study of acoustic oscillations in the Sun - allows us to probe solar interior structure and rotation in 
two dimensions, depth and latitude, by taking advantage of the different spatial distribution of the various modes. Recently, 
two projects - the Global Oscillation Network Group (GONG) and the Solar Oscillations Investigation using the Michelson 
Doppler Imager (MDI) instrument aboard the SOHO spacecraft - have provided nearly five years of continuous helioseismic 
data, allowing the rise of the current solar cycle to be followed in unprecedented detail. 

Two-dimensional inversions for the rotation profile have been widely studied. Inversions for the structural parameters (for 
example, sound speed and density) are more commonly carried out in one dimension only. However, it is possible to access 
two-dimensional information on these parameters also, using the so-called 'even a coefficients'; results of early attempts can be 
found in Gough et al. (1996). In the light of recent work by Howe et al. (2000h), Howe, Komm and Hill ( ^OOOa ) , Toomre et al. 
(2000) and Antia and Basu (200C) suggesting that the torsional-oscillation pattern seen at the surface penetrates substantially 



into the convection zone, it is of interest to see whether the structural changes also penetrate deeper than has previously been 
thought. In the work presented here we examine changes in helioseismic determinations of the solar asphericity as the solar 
cycle progresses. 

The modes are described by the radial order n, related to the number of nodes in the radial direction, and the degree I 
and azimuthal order m which characterize the spherical harmonic defining the horizontal structure of the mode. Rotation and 
asphericity lift the degeneracy between the 2l-\-l modes of different m making up an (n, I, m) multiplet, resulting in frequency 
splitting. The frequencies Unim of the modes within a multiplet can be expressed as an expansion in orthogonal polynomials, 
for example 



(1) 



In early work (e.g., Duvall, Harvey & Pomerantz, 1986) Legendre polynomials were commonly used, whereas more recent 
work often uses the RitzwoUer-Lavely formulation of the Clebsch-Gordan expansion, where the basis functions are polynomials 



related to the Clebsch-Gordan coefficients (RitzwoUer fc Lavely 1991) CjJJ|„ by 
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V(2/-j)!(2/ + j + l)! 



are referred to as a-coefEcients. The odd-order coefficients, 01,03,..., are used to 
determine the rate of rotation inside the Sun, and reflect the advective, latitudinally symmetric part of the perturbations 
caused by rotation. The even coefficients, which are much smaller, are sensitive to second order contributions from rotation, 
any possible magnetic field and any possible departure of the solar structure from the spherically symmetric state. Since the 
rotation rate in the solar interior is determined by the odd-order coefficients, it can be used to calculate the second order 
contribution of rotation to the even-order coefficients. Further, since the time variation in the rotation rate is rather small, 
being of the order of 0.5%, the second order contribution due to rotation may be expected to be essentially constant in time 
and hence any time- variation in the even-order coefficients should be due to other sources, e.g., magnetic field or asphericity. 
Unfortunately, it is not possible to distinguish between these two possibilities using the even-order coefficients (Zweibel and 
Gough 1995). 

It was well established during the previous solar cycle, that the even coefficients show temporal variations related to 
solar activity measures. The 02 and 04 coefficients of Duvall et al. (1986) were used by Gough and Thompson (1988) to infer 
the possibility of a shallow magnetic perturbation in the sound speed near the equator. Kuhn (1988a) also pointed out the 
relationship between 02 and 04 measurements from BBSO and earlier data and the latitudinal variation of internal sound 
speed, noting the possibility of temporal variation as the known 'hot bands' associated with magnetic activity migrated during 
the solar cycle. The temporal variation of the 02 and 04 coefficients and its relation to changes in the latitudinal dependence of 
limb temperature measurements was further studied by Kuhn (1988b), Libbrecht (1989) and by Kuhn (1989), who predicted 
that the same relationship should extend to higher-order coefficients. The inference from the inversions of the BBSO even 
a coefficient data carried out by Libbrecht and Woodard (1990) and Woodard and Libbrecht (1993) was that most of the 
variation in the even coefficients was localized close to the surface and at the active latitudes, with a near-polar variation 
anticorrelated to the global activity level. In the new cycle, we have data from GONG and MDI which allow us to study 
these trends in more detail. Howe, Komm and Hill (1999), have found linear relations between the even-order a coefficients 
and the Legendre decomposition components of the surface magnetic flux up to as, and showed that these relations could be 
extended backwards in time to the BBSO data from the previous cycle. Dziembowski et al. (2000) found a good correlation 
between the even-a coefficients derived from 12 72-day MDI data sets (1 May 1996 - 31 May 1998) and the corresponding 
even components of the Ca II K data from Big Bear Solar Observatory (BBSO) up to aio. 

It has been pointed out by Kuhn (1998) that the variations in mode frequency and even a coefficients are unlikely to 
be caused directly by the surface flelds, as the flelds required would be substantially stronger than those observed. However, 
this does not rule out the idea that the presence of magnetic flelds affects the temperature. As pointed out by Dziembowski 
et al. (2000), the magnetic held might influence the thermal structure through the annihilation of the fleld (/3 effect) and 
indirectly through the perturbation of the convective transport (a effect) in addition to the direct mechanical effect of the 
Lorentz force. The good correlation between the magnetic activity indices and the variation in mean frequencies and even 
a-coefficients would tend to suggest that magnetic fleld probably plays some role in these temporal variations. 

In the present work we use more extensive data sets from both GONG and MDI, extending the GONG analysis up 
to ai4, and present the results of 2-dimensional inversions for the sound speed. We effectively assume that the even-order 
a-coefficients arise from aspherical sound speed distribution, rather than being due to possible magnetic fleld effects. 



2 DATA 

We have analysed 53 overlapping 108-day time series of GONG data covering the period 1995 May 7 to 2000 October 6, 



centered on dates 1 36-day 'GONG month' apart. The data were analysed through the standard GONG pipeline (Hill et al 



199f)), to flnd frequencies for each mode, and then a-coefficients were derived by fltting to the frequencies in each {n,l) 
multiplet. A typical mode set from this process would consist of around 1400 multiplets with coefficients up to ais, for 
I < 150; the intersection of all the sets contains about 400 multiplets mostly with I < 100. 

The MDI data consist of 23 non-overlapping time series covering the period 1996 May 1 to 2001 April 4, although with 



some interruptions due to problems with the SOHO satellite. The data were analysed as described by Schou (1999). These 
sets contain coefficients, up to aae, for roughly 1800 multiplets with I < 300. The substantial temporal overlap between the 
two sets allows useful cross-checking of results. 

For the purpose of the inversions we have used only p-mode splittings for modes with frequencies between 1.5 to 3.5 mHz. 
Further, only coefficients up to ai4 are used as the higher order coefficients do not appear to be signiflcant. This gives us 
typically 8000 coefficients for MDI data sets and 6000 coefficients for the GONG data sets. Frequencies for f-modes {n = 0) 
are also available for the MDI data, and might provide additional information, but they have not been used in the current 
work. 
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3 COEFFICIENT ANALYSIS 



3.1 Method 



Howe et al. ( 1999 ) considered the temporal variation of the GONG mean even a coefficients up to as and showed that they 
were strongly correlated with the corresponding components of the Legendre decomposition of the surface magnetic field. In 
the present work we extend this analysis, with some modifications, to coefficients up to ai4 and to MDl as well as GONG 
data. 

The latitudinal distribution of the surface magnetic fiux can be expressed as a sum of Legendre polynomials in cos d, 



B{t, cosd) ^^Bk{t)Pk{cose). 



(3) 



where B^it) axe time- varying coefficients. These coefficients can be compared with the corresponding frequency splitting 
coefficients. Global helioseismic measurements are sensitive only to the latitudinally symmetric part of any departures from 
spherical symmetry — that is, to the even components of the expansion. 

For the purposes of the present work, we consider two weighting schemes. In the first, the splitting coefficients were 



weighted by Eni, where Eni is the mode inertia defined as in Howe et al. (1999) and references therein, normalized to set the 
value for the radial mode at 3 mHz to unity. In the second they were weighted by lEni/Qik, where Qik is defined in eq. pT| 
below, before averaging over all common multiplets. This removes the main /-dependent part of the coefficient variation, and 
also reduces the frequency dependent part. Thus, we define the quantities 



< a2k > (i) = X! 



a2fc(n, /, t)E„i 



n,l 



and 



< hk >{t) = J2 



a2k{n,l,t)lE„i/Qik 



PEl,al,, (n, Z, t)/Q^^^ ' ^ PEl.al,, (n, Z, t) ' 



(4) 



(5) 



where t is time and Oa^ is the estimated uncertainty in the coefficient ak, and the sum is over the approximately 400 {n,l) 
multiplets that are common to all the mode sets. 

We then express the variation of the < ak > and < 6^ > as a function of the Bk, for even k, as 



< Ofe > (t) = Cfc ~f mkBk{t) 



(6) 



and perform linear least-squares fits to obtain the gradient nik and intercept Ck for each k. The intercept will contain a 
contribution from the invariant part of the coefficients which includes the second order effect of rotation, as pointed out by 



Woodard (1989), Gough and Thompson (199C), and Dziembowski and Goode (199S 



3.2 Results 

In Figure |l| we show the variations of the mean coefficients < 02 > to < 014 > as a function of time, for the sets of modes 
common to all datasets. The correlation between the coefficients and the Legendre components of the surface fiux persists up 
to ai4 in both sets of data (Table |l|). (The value of the correlation coefficient R required for the 0.1 per cent significance level 
is 0.45 for GONG and 0.68 for MDI.) There is a reasonably good agreement between GONG and MDI data for 02, . . . ,112, 
but for ai4 the agreement is less good. We suspect the difference may be due to the different analysis procedures used for 
calculating the splitting coefficients in the two data sets. In principle. Bo may be correlated to the temporal variations in the 
mean frequencies, but in practice, it is difficult to separate out the contributions to the mean frequency due to the magnetic 
field from those due to uncertainties in the spherically symmetric structure of the Sun. For this reason, we do not compare 
the component Bo in this work. 

In Figure ^ we show the coefficients < 02 > to < ai4 > as a function of the corresponding components B2 to B14 of the 
surface magnetic flux, for GONG and MDI, while the left portion of Table |l] gives the slopes and correlation coefficients 
for the fits between even a-coefficients and their corresponding magnetic flux decomposition components. While alternate 
even coefficients < 02 >,< aa >,< oio >, . . . appear to be anticorrelated with Bk, the other coefficients < 04 >, < ag >, . . . 
are correlated with Bk- There also appears to be a reduction in magnitude of < > with k for k > 10, while the Bk do 
not show such a marked decrease. However, the < b2k >, (Figure ^, right side of Table ^ are all positively correlated with 
B2fc with scaling constants that appear to be independent of k, except for the highest order coefficients where the GONG 
data show some reduction in sensitivity. In fact, most of the variation, including the sign changes, in the < ak > slopes up to 
k = 12 can be explained by the corresponding variation in the angular integrals, Qik- 
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Figure 1. Temporal variation of the GONG (open circles) and MDI (filled triangles) mean even a coefficients. The quantity plotted is 
< a^. >, as defined in Eq. ^. The curves show the best-fit values [m^Bf. + Cf.] for linear fits between the < > coefficients and the 
corresponding B2k components of the Legendre decomposition of the magnetic flux, for GONG (solid) and MDI (dashed). 
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Figure 2. The relations between the GONG (open circles) and MDI (filled triangles) < aj. >-coeflicients and the corresponding Bf^ 
components of the Legendre decomposition of the magnetic flux from the Kitt Peak synoptic maps. The lines show the best-flt results 
for linecir fits to the data, for GONG (solid) and MDI (dashed). 
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Figure 3. The relations between the GONG (open circles) and MDI (filled triangles) < &j. >-coefiicients and the corresponding B^. 
components of the Legendre decomposition of the magnetic flux from the Kitt PeaJj synoptic maps. The lines show the best-fit results 
for lineaj: fits to the data, for GONG (solid) and MDI (dashed). 
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4 ASPHERICITY INVERSIONS 
4.1 Inversion techniques 

In order to study the variation of aspiiericity witli deptii and latitude we apply an inversion technique to the even-order 
splitting coefficients. We use the variational principle to analyse the departures from a spherically symmetric solar model 



(e.g., Gough (199J)), in order to study aspherical perturbations to the sound speed and density in the solar interior. For 
simplicity, we only consider axisymmetric perturbations (with the symmetry axis coinciding with the rotation axis) that are 
symmetric about the equator. In this case, using the variational principle, the difference in frequency between the Sun and a 
solar model for a mode of a given order, degree and azimuthal order (n, /, and m) can be written as: 

^ = r dr r r sinOde ^(^r)^^{r,e)+fC;'Ar)^-{r,e)\ YriYrT (7) 

fulm Jo Jo Jo \ P J 

where r is radius, Q is colatitude, 5v„im / fnhn is the relative frequency difference, IC^l p{r) and /C"'j,2(r) are the kernels for 



spherically symmetric perturbations (Antia and Basu 1994), and Y/" are spherical harmonics denoting the angular dependence 



of the eigenfunctions for a spherically symmetric star. It is assumed that the Yj™ 's are normalized such that 

d(i> / sine de YriYry ^ I. (s) 

Jo 

The perturbations 5<? ji? and 5p/p can be expanded in terms of even order Legendre polynomials: 

^(r,e)=^Cfe(r)P2fe(cose), (9) 

^(r,0) = y pfc(r)P2fe(cos6)), (10) 

where Cfc(r) and Ph(f) are shorthand notations for (i5c'^/c'^)i,(r) and {5p/ p)k{r) respectively. The spherically symmetric com- 
ponent (fc = 0) gives frequency differences that are independent of m and thus only contribute to the mean frequency of the 
{n, I) multiplet. Higher order terms give frequencies that are functions of m and thus contribute to the splitting coefficients. 
The angular integrals in Eq. ^ can be evaluated to give 



d(t> / sm9deYr{YrYP2k{cose) = -QikV},>{m) (11) 



where Qik depends only on I, k and 'P2fc ("^) orthogonal polynomials defined by Eq. |^. The extra factor of l/l ensures 

that Qik approach a constant value at large Thus with this choice of expansion (Eq. ^) the inversion problem is decomposed 
into independent inversions for each even splitting coefficient and Ck{r), pk{r) can be computed by inverting the splitting 
coefficient a2k- This would be similar to the 1.5d inversion to determine the rotation rate (RitzwoUer & Lavely 1991), called 
so because a two dimensional solution is obtained as a series of one dimensional inversions. 

In practice, we also need to account for the contribution to the frequency splittings that arises from uncertainties in 
the treatment of surface layers in the model. It is known that this contribution to frequency splittings should be a slowly 
varying function of frequency alone once it is corrected for differences in the mode mass of the modes; no degree dependence 
is expected to a first approximation. As in the case of inversions to determine the spherically symmetric structure of the Sun, 
the surface uncertainties are accounted for by assuming that these can be represented by a frequency-dependent function 
F{u). However, in this case we determine a different function Fkiu) for each coefficient a2k and write the inversion problem 
as 

la2k[n,t) ^ / Fk{l^, 
= Qik / JC^2 „Ck(r) dr + Qik / IC 2pk{r) dr + Qik — — 

Vr,l Jo Jo Fr, 

Eq. |l^ can be inverted using the usual inversion techniques to calculate Ck{r), pk{r) and Fk{y). We use a Regularized Least 
Squares (RLS) technique for this purpose and expand the unknown functions Ck (r) , pk (r) , Fk (v) in terms of cubic B-spline basis 
functions over knots, which are approximately uniformly spaced in acoustic depth (or frequency). First derivative smoothing 
is applied to constrain the error in the solution to remain small. 

As in the case of inverting for rotation rate it is possible to apply a 2d inversion technique for asphericity, where the 
asphericity is not expanded in terms of Legendre polynomials, but instead one uses 

la2k{n,l) Qik{4.k + 1) f'^ H . „ ,„ f y.„i 5c^ Sp\ /)^ , ^ Fk{iy„t) , , 

- I dr smO de [ K:^2 „^ + IC ^2— ] P2k{cose) + Qik — - — - (13) 



JC:ipCk{r) dr + Qik / JC;'2Pk{r) dr + Qik^^ (12) 
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Table 1 . Gradients m and correlation coefBcients R for fits between the weighted mean coefficients <aj.>,<fefe> and the corresponding 
Legendre components of the unsigned magnetic flux. 





GONG Ofc 




MDl CLk 




GONG bk 




MDI bk 




k 


rrij. (nHz) 


R 


rrik (nHz) 


R^ 


rrik (nHz) 


R 


rrik (nHz) 


R 


2 


-0.18 ±0.00 


-1.00 


-0.18 ±0.00 


-1.00 


24.81 ±0.13 


1.00 


25.76 ±0.16 


1.00 


4 


0.14 ±0.00 


0.98 


0.13 ±0.00 


0.96 


25.57 ±0.36 


0.98 


26.50 ±0.57 


0.96 


6 


-0.11 ±0.00 


-0.99 


-0.10 ±0.00 


-0.99 


26.50 ±0.26 


0.99 


26.45 ± 0.43 


0.99 


8 


0.09 ± 0.00 


0.99 


0.08 ±0.00 


0.99 


27.06 ±0.33 


0.99 


27.00 ±0.51 


0.99 


fO 


-0.06 ±0.00 


-0.99 


-0.06 ±0.00 


-0.99 


25.70 ±0.46 


0.99 


25.83 ±0.64 


0.99 


12 


0.05 ±0.00 


0.97 


0.05 ±0.00 


0.98 


24.40 ±0.79 


0.97 


24.81 ±0.98 


0.98 


14 


-0.02 ±0.00 


-0.93 


-0.04 ±0.00 


-0.97 


14.56 ±0.64 


0.93 


27.67 ± 1.30 


0.97 



where &c? jc? and (5p/p are now functions of (r, and can be expanded in terms of a set of 2d basis functions for the RLS 
inversion. The knots are chosen to be uniformly spaced in acoustic depth and cos(S) and we use the product of B-spline basis 
functions in r and cos(6'). Thus, for example, we can write 



„2 



where hij are the coefficients of expansion and <i>i(r) are the B-spline basis functions over r and ^j(cos9) are those over 
cos 6; Ur and ng axe the number of basis functions in r and cos 6, respectively. The unknown functions — and Fk{u) 
axe expanded in terms of basis functions and the coefficients of expansion are determined by solving one inversion problem 
involving all splitting coefficients. First derivative smoothing is applied in both r and cos 6. Thus, the coefficients of expansion 
are determined by minimizing 

, o r , ^ ,x /. R /.I / ,o 

-I ■ / / acosc I /v,2 „ir,Pj — ^r,(:?j ,2ir,pj-^ 



E 



QlkFk(v) 
Enl 



+ A6 



R p1 

dr d cos 
J -1 



d_ I 5c- 

dr 



dr I dcosS sin 



d 



9 cos 6 V 



+ 



5c' 



d_ fsp 

dr \ p 



d 



Sp 

d cos 6 \ p 



(15) 



where Ar and Xg are the two regularization parameters controlling the smoothing and and /C" *2 are the 2d kernels, (t„ 

are the uncertainties in a2k{n,l). We have used 16 knots in r and 10 knots in cos 6 to represent the asphericity. 



4.2 Inversion Results 

In Figure ^ we show the results of sound speed inversions of GONG and MDI data at selected radii as a function of time and 
latitude. There is no signature of significant temporal evolution at these radii. Both sets of data show a persistent sound-speed 
excess at around 60 degrees, apparently extending well down into the convection zone, though this is more significant in the 
MDI than in the GONG data. 

The inversions show the deviation from a spherically symmetric model. It turns out that the temporal mean of these 
deviations is non-zero and has a dependence on depth and latitude, as illustrated in Figure ^. However, it may be noted that 
the current data extend over only a small fraction of a magnetic cycle and the temporal mean over this limited period may 
not have much significance. The second order contribution from rotation would also contribute to the temporal mean. There 
is some difference between the mean calculated from GONG and MDI data sets, which could be due to systematic differences 
between the two sets. In both cases there is a peak around depth of 0.08-Rq and the peak is more pronounced in MDI data. 

The sound-speed variation after subtraction of the mean profile is shown in Figure ^ No systematic structure is evident. 
Figure ^ shows the variation of the sound-speed residuals from the 2d inversions of GONG and MDI data at selected (r, 6) 
points, illustrating the agreement between the two experiments. The residuals are obtained by subtracting the temporal mean 
from the results at each epoch. It is clear from these figures that there is no significant temporal variation in the asphericity 
at the depths we can resolve. 

The strong variation in the even a coefficients is evidently not reflected in the inversion results. Instead, it has been 
absorbed in the surface terms. To illustrate this, we show in Figure ^ the reconstruction of the surface term for a frequency 
of 2.5 mHz, with overlaid contours of the magnetic flux. As we would expect, since the dependence of the individual surface 
terms on frequency and on the magnetic field strength is largely independent of the order of the coefficient, this reconstruction 
matches the magnetic flux quite well. This finding is consistent with the results of Woodard and Libbrecht (1993), who carried 
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Figure 4. Grey scale maps showing the results of sound-speed inversions of GONG (first and third columns) and MDI (second and 
fourth columns) data at radii 0.98R (top), 0.88il (centre) and O.SOi? (bottom). Blank spaces represent periods where there are no data 
available. Columns 1 and 2 show the 1.5d inversion results and columns 3 and 4 the 2d results. 
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Figure 5. Grey scale maps showing the mean results of sound-speed inversions of GONG (left) and MDI (right) data, multiplied by 
10*, as a function of latitude and radius. 



out latitudinal inversions of the coefficients up to 012. They found that the sound speed showed a peak at the active latitudes, 

matching both magnetic and temperature behaviour, and also a decline at high latitudes with increasing solar activity which 
was better matched by the temperature than by the magnetic data. The surface terms from our inversion show a similar high- 
latitude behaviour, which does in fact correspond well to the high-latitude magnetic data; we attribute this to an improvement 
in the magnetic observations since 1992. Our inference that most of the solar-cycle variation in the sound speed is localized 
near the surface is also in agreement with the conclusions of Woodard and Libbrecht (1993). 
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Figure 6. Grey scale maps showing the results of sound-speed inversions of GONG (first and third columns) and MDI (second and 
fourth columns) data at radii 0.98/? (top), 0.88-R (centre) and 0.80-R (bottom), after subtraction of the temporal mean. Blank spaces 
represent periods where there are no data available. Columns 1 and 2 show the 1.5d inversion results and columns 3 and 4 the 2d results. 



5 RESOLUTION ISSUES 

In order to test our inversion method and to see how well the inversions can resolve different features, we have conducted tests 
with artificial data. Since asphericity inversions do not include the spherically symmetric component corresponding to the 
mean frequencies, &<?!<? integrated over a surface with constant r must vanish. This arises because the non- vanishing part of 
the integral can only contribute to the mean frequencies which are not included. Furthermore, the even splitting coefficients 
are sensitive only to the north-south symmetric component of the asphericity. Thus it is necessary to choose artificial data 
which assume a profile that integrates to zero in latitude and which are symmetric about the equator. We have chosen a 
profile of the form 

|aoexp[-(^)2]p3(cos2 0) if jr - r„l < d (^g^ 
c otherwise 

where ro and d are constants which respectively define the mean depth and half-width of the aspherical perturbation and 
P3(cos^ ff) is a third degree polynomial in cos^ Q which integrates to zero. There is no special reason to choose a third degree 
polynomial, except that with this choice one can get two bands of positive and negative asphericity in each hemisphere. The 
amplitude ao of the signal is chosen to approximately match the observed signal. 

Using the artificial profile for 5(? /(? with the peak located at different radii and with different widths, we construct 
a set of artificial data. These data include only those modes which are present in the observed data, and random errors 
with a Gaussian distribution with standard deviation equal to the estimated uncertainties in observed data are added to the 
calculated splittings. The amplitude ao is chosen to be comparable to that obtained from real solar data. This process is 
followed for both GONG and MDI data sets. These artificial data were then inverted the way we would invert the real data 
and using the same choice for smoothing as is used for real data. The inversion results can be compared with the actual profile 
used in constructing the artificial data. 

The inversion results show that we are able to invert for the aspherical sound-speed difference in the outer layers of the 
Sun. Fig. ^ shows the inversion results for artificial data with peak located at different depths. These results indicate that 
with the current data sets we cannot invert for any feature located below a radius of 0.8 Rq. At smaller radii, e.g., Q.IRq, 
although we get some result at the equator, the artificial profile is not reproduced properly at high latitudes. Similarly, the 
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Figure 7. Variation of GONG (open) and MDI(filled) sound-speed residuals at radii 0.98ii (top), 0.86R (middle), and 0.8R (bottom) 
for latitudes 0° (left) and 20° (right.) 
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Figure 8. Grey scale map showing the reconstruction of the latitudinal dependence of the surface term from GONG (left) and MDI 
(right), multiplied by 100. Overlaid contours show the Kitt Peak unsigned magnetic flux with the Bo term subtracted; contour spacing 
is lOG. 



results become unreliable above 0.98 Rq, mainly because of the lack of high-degree modes. The inversion results are reliable 
between depths of 0.027?© and about O.2J?0. As with all RLS inversions, there is some structure in the results away from the 
peak, but the structure is small and does not give rise to any major artifacts. The figures shown give us confidence in the 
results we obtain by inverting real data. 

A similar exercise was also carried out for the inversions of density, but with an amplitude of artificial profile comparable 
to those in real data, it was not found to be feasible to reproduce the profile to any reasonable accuracy at all depths. Thus 
it appears that the aspherical component of density cannot be reliably determined from current data sets. Therefore, we do 
not show these results in the present work. 
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Figure 9. Artificial data inversion results for a variety of radial distances (from top to bottom, ro = 0.7R, 0.8R, 0.9R, 0.95R, 0.98i?) 
and latitudes (from left to right, latitude 0°, 30° and 56°), for GONG (dotted curves) and MDI (dashed curves.) The solid curves show 
the 'true' rate and the thinner dashed and dotted curves represent l-a error bounds on the inversion results. 



6 DISCUSSION 

We have studied the temporal variations in even a-coefiicients from GONG and MDI data during the period 1995-2001. The 
mean values of the even o-coefHcients (after scaling for mode mass and angular integrals) over all modes are found to be well 
correlated to the corresponding components of the observed magnetic flux at the solar surface. The slope of best linear fit 
between the mean splitting coefficient and the corresponding Legendre component of the surface magnetic flux is found to be 
essentially independent of the order of the coefficient. Thus the latitudinal variation in the surface magnetic flux is correlated 
with that in asphericity as measured by the even a-cocfficicnts. Dziombowski ct al. (2000) obtained similar results analyzing 
12 72-day MDI data sets (1 May 1996 - 31 May 1998). They found a good correlation between the even-a coefficients (up to 
aio) and the corresponding even components of the Ca II K data from Big Bear Solar Observatory (BBSO). 

In this work we have done two dimensional inversions for sound speed using even arcoefficients from GONG and MDI 
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data sets covering the period from 1995 to early 2001, encompassing the rising phase of the current solar cycle. We find 
no significant temporal variation in the asphericity of the sound speed over this period. Thus it appears that the temporal 
variations in even a-coefficients arise from changes taking place in surface layers. A similar conclusion has been obtained for 
the spherically symmetric component of sound speed and density (Basu and Antia 2000) as the changes in mean frequency 
also appear to be associated with surface effects. 

For this work we have assumed that the even a-coefficients arise from aspherical sound-speed distribution, but this is 
by no means obvious as these coefficients could also arise from a magnetic field since both sets of kernels are very similar in 
most parts of the Sun and cannot be easily distinguished (see, for example. Figure 10 in Dziembowski et al. 2000). In that 
case a magnetic field strength yielding v\/c? ~ 10~* would be required, where va is the Alfven speed. Inside the convection 
zone one might not expect an ordered magnetic field over large length scales as turbulence might be expected to randomize 
the magnetic field. Such a randomized magnetic field can also effectively change the wave propagation speed, giving a signal 
similar to that from aspherical sound-speed distribution. It may not be possible to distinguish between these two effects from 
the even a-coefficients. Kuhn (1998) has argued that the observed magnetic field at the solar surface is not sufficient to explain 
the magnitude of the even a-coefficients. However, one can argue that the magnetic field increases significantly as one goes 
deeper in the near surface layers and this could explain the observed magnitudes of the splittings. 

The aspherical component of the density cannot be reliably determined from the current data sets, but the magnitude is 
generally found to be smaller than that for the sound speed. If the aspherical perturbations were of thermal origin then one 
might have expected the density perturbations to be comparable to the sound-speed perturbations. 

Although the temporal variation in asphericity in the solar interior is not significant, the mean asphericity is found to 
be significant and shows a peak around r = O.92_R0 which is very pronounced in the inversions of MDI data. This is deeper 
than the depth to which the surface shear layer in solar rotation rate extends (Schou et al. 1998) and about the same as the 
penetration depth recently (Howe et al. 2000; Antia and Basu 2000) established for the torsional oscillation pattern. Antia, 
Chitre and Thompson (2000) attempted to detect a possible signature for magnetic field using the even a-coefficients to 
find that the coefficients a2 and 04 have a residual after subtracting the contribution from rotation, which is peaked around 
r — 0.95Rq, which is consistent with our results. Dziembowski et al. (2000) calculated inversions separately for temperature 
and magnetic field perturbations and detected a significant perturbation in the spherically symmetric part at a depth of 
25-100 Mm with a maximum at about 45 Mm, which agrees with the depth at which we find the aspherical perturbation 
in the sound speed. They interpret this perturbation as being due to either a magnetic perturbation of about (60 kG)^ or a 
temperature perturbation of about 1.2 10"'* which is of the same order as the result presented here. Oilman (2000) pointed 
out that a thermal driving mechanism for the observed meridional fiow would require a temperature excess at the equator — 
the opposite of what we observe. Thus this finding poses yet another challenge for theoretical understanding of the dynamics 
of the convection zone. 
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